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Abstract 



In order to understand the long known anomalies in the composition dependence 
of diffusion and viscosity of binary mixtures, we introduce here two new models 
and carry out extensive molecular dynamics simulations. In these models, the two 
molecular species (A and B) have the same diameter and mass. In model I the 
inter-species interaction is more attractive than that between the pure components, 
while the reverse is true for model II. Simulations and also mode coupling theory 
calculations reveal that the models can capture a wide variety of behavior observed 
in experiments, most interesting among them are the non-monotonic variation of 
diffusion and viscosity with the composition and the re-entrant viscosity dependence 
of the relaxation time. 
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Transport properties of binary mixtures often show strong and baffling depen- 
dencies on the composition which have not been understood or even adequately 
addressed to in a molecular theory. The well-known Raoult's law of classical phys- 
ical chemistry^ predicts the following simple hnear dependence on the composition 
for a given property P, 

P = XiPi + X2P2 (1) 

where the mole fractions and PjS are the values of the property P of the 

pure (single component) liquids. More often than not, significant deviation from 
Eq. 1 is observed. Of many anomalies exhibited by binary mixtures, the existence 
of an extremum (sometimes even a double extremal) in the composition dependence 
of excess viscosity^ and the re-entrant type behavior of the relaxation time when 
plotted against viscosity^ are certainly the most remarkable. The latter shows, in 
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a dramatic fashion, that viscosity is not a unique determinant of relaxation in bi- 
nary mixtures'^. Several interesting theoretical and computer simulation studies on 
Lennard- Jones binary mixtures have been carried out recently but these studies 
have concentrated mainly on the glass transition in binary mixtures and considered 
only one particular composition and a unique interaction strength. Earlier Heyes 
carried out the extensive equilibrium MD simulations of Lennard- Jones binary mix- 
tures by using both the microcanonical (N V E) and canonical (N V T) ensemble 
methods to study the partial properties of transport coefficients in the inert gas 
medium ^. The non-ideality in the case of inert gas mixtures is small, since their 
mutual interaction strength {eAs) follows the Berthelot mixing rule. 

In order to understand the markedly nonlinear composition dependence, here we 
introduce and study two models (referred to as model I and model II) of binary 
mixtures in which the solute-solvent interaction strength is varied by keeping all the 
other parameters unchanged. In our models, all the three interactions (solute-solute, 
solvent-solvent and solute-solvent ) are described by the (6 — 12) Lennard- Jones 
potential. 



where i and j denote any two different particles. We set the diameter (a) and 
mass (m) of both the solute and the solvent molecule to unity, for simplicity. The 
solvent-solute interaction strength lies in the potential well depth eAB, where A and 
B represent the solvent and solute particles, respectively. Throughout this study we 
keep the interaction strength eaa = 1-0, (solvent-solvent), ess = 0.5 (solute-solute). 
In the two models we use two different solvent-solute interaction strength values, 
namely cab — 2.0 in model I and cab — 0.3 in model II. So, while model I is a 
"structure former" (between A and B), model II is a "structure breaker". Note 
that both A and B have the same diameter. We believe that these simple models 
can serve as starting points to understand the many baffling properties of binary 
mixtures. 

Extensive MD (microcanonical ensemble, with the usual periodic boundary con- 
dition) simulations have been carried out with a total of 500 particles for two types 
of models by varying the solute mole fraction (that is, of B) from to 1. The re- 
duced temperature T* {— kBT/e) is set equal to unity in model I and 1.24 in model 
II and the reduced density (p* = pa^) is 0.85 in both the models. After many trial 
runs to verify the existing results on viscosity of one component liquids, we have 
selected a time step At* = 0.002r for model I and At* = O.OOlr for model II. where, 
r = a^m/e. We have dealt with six different solute compositions, namely 0, 0.2, 
0.4, 0.6, 0.8 and 1.0. For each solute composition we have equilibrated the system 
up to 1.5 X 10^ steps. Simulations carried out for another 2.0 x 10^ steps after the 
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equilibration during which all the relevant quantities have been calculated. For each 
composition, we have run three independent simulations and have taken averages 
over them. We have checked all the three partial radial distribution functions to 
make sure no phase separation occurs during simulations (for model II). Viscosity 
values are obtained by integrating the stress time correlation function which defines 
the time dependent viscosity by the following relation 

r){t) = {VkBT)-\a^^{0)a^^{t)) (3) 



where the off-diagonal element of the stress tensor a^^ for binary mixture is defined 
as, 

Ni N 

<^"^ = E[bM/H + ^^^.]+ E UpMM + ^I^^] (4) 

i=i j=Ni+i 

Here, Fj is the z-component of the force acting on the j-th particle and the cor- 
responding X-coordinate is xj, pj is the z-component of the momentum pj of j-th 
particle, m being the mass of the particle. Amongst total N particles, Ni solvent 
particles are labelled from 1 to A^i and solute particles from {Ni + 1) to N. 

Diffusion coefficients are calculated both from the mean square displacement and 
velocity autocorrelation function via the Grccn-Kubo relation. The results of the 
simulation are given in figures 1 — 4. We shall discuss the results after we describe 
the mode coupling theory employed. 

We have carried out mode coupling theory (MCT) calculations of diffusion 
and viscosity to understand the simulation results, especially the origin of non- 
monotonicity. These calculations have been carried out by using well-estabfished 
expressions ^"^^ ]\fQ^g that for binary mixtures at normal density and temperatiu^c, 
the short time dynamics of the relevant time correlation functions are important 
and in fact, can contribute more than 50% of the total value, just as for one com- 
ponent liquids. Thus, any solution of the MCT equations require accurate input of 
the short time dynamics. For a given transport property P, MCT formalism^"^^ 
assumes the following separation into the short time, binary collision controlled, 
contribution p(^'") (t) and the contribution from the collective term, which in dense 
fiquid is dominated by the density term, P'^PP\t). So the total dynamical quantity 
P{t) can be written as ^, 

P{t) = P^^^\t) + P^pp\t) (5) 
The binary part of both viscosity and the friction are assumed to be Gaussian, 

p(''-)(t) = p{t = 0)exp{-ty4). (6) 
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For viscosity, P{t — 0) is the high frequency shear modulus Goo of a binary mixture 
given by, 



Goo = {pi + p2)kBT + ^ X! PiPj J ^^9ij{r)-^ 



dr 



(7) 



here, i,j = l indicate solvent particles and i,j = 2 denote solute particles, pi and 
P2 are the coresponding number densities for the solvent and solute particles. gij{r) 
is the partial radial distribution function of the particles labelled i and j. In the 
present calculation, gij{r) is obtained from Ornstein-Zernike equations with SMSA 
closure^^, which provides a reasonable, although not perfect, agreement with the 
simulation results. 

The characteristic relaxation time r^j for viscosity can be determined by the 
second derivative of ri[t). The resulting expression of for a binary mixture is rather 
complex and is given elsewhere The mode coupling contribution to viscosity is 
assumed to be given by the binary product of the density terms^'^. In the present 
case, one derives contribution from four such density terms and total MCT term 
jy(pp) is given by the simple addition of four rj'^PiPi') terms. The frequency dependent 
diffusion coefficients Di{z) are related to the respective frequency dependent frictions 
according to Einstein relation. 

For friction, Q[t = 0) is the Einstein frequency of the i-th component in the mixture 
and is determined by the static correlation functions. The initial short time part 
of time dependent friction is assumed to be Gaussian with characteristic time 
which is calculated using the short time expansion of the force-force time correlation 
function The expression for the mode coupling term has already been given 

by Bosse et al and need not be presented here. The dynamical input parameters 
for MCT calculation are the wavenumber (q) dependent and time (t) dependent 
partial intermediate scattering functions, Fij{q,t) and the self-dynamic structure 
factors, Fsi{q, t). The expressions of dynamic structure factors Fij{q, z) are obtained 
by solving four coupled equations obtained from time dependent density functional 
theory This method requires the value of the frequency dependent self-diffusion 
coefficient of both the species. We have used the bare value of the diffusion coefficient 
here, calculated from Eq. 8 by using the binary part of the frequency dependent 
friction. We approximate Fsi{q,t) by Fsi{q,t) — exp{ ~'^ where 

{Ar^{t)) = 2 f C,i{T){t - T)dT (9) 
Cyi{t) in Eq. 9 is obtained from the numerical inverse Laplace transform of Cv{z). 
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Figure 1 shows the remarkable re-entrant behavior of the structural relaxation 
times (tj) when the viscosity is changed by varying the composition. The increase in 
composition is indicated by arrows. Here the relaxation time Tj is calculated by using 
Tj = (T^/Dj. Thus, the relaxation time is inversely proportional to diffusion constant. 
We believe that a similar behavior will be observed for rotational relaxation as well. 
The simulation points here are averages over three independent long runs; error bars 
are typically ±0.2 for viscosity and ±1.5 for the relaxation time (that is, ±0.002 for 
diffusion coefficients). Note that in this figure we have shown only the simulation 
results, for clarity - theory shows a similar behavior. 

Figure 2 shows the composition dependence of diffusion coefficients obtained 
from both theory and simulation, for model I. Figure 3 shows the same for model 
II. Note the non-monotoic composition dependence. Diffusion of A and B show 
differing behavior, in all the cases. Figure 4 depicts the nonideality of viscosity 
with respect to composition, for both the models. Though the agreement between 
theory and simulation is certainly not perfect, the trends are similar in both the 
calculations. Note that the theoretical calculation does not use any simulation data 
as input or any adjustable parameter either; thus the theory and the simulation 
provide independent test of each other which is important for binary mixtures. 

We conclude this Letter with the following comments. 

(1) It is shown that the simple models can describe the decoupling of diffusion 
from viscosity in binary mixtures, when the viscosity is changed by varying the 
composition. This decoupling is most dramatically manifested in the re-entrant type 
behavior depicted in Fig. 1. This shows that viscosity is not a unique deteminant 
of relaxation in binary mixtures. 

(2) Models I and II seem to reproduce the behavior observed in large number 
of systems. We believe that this is the first time a microscopic model captures the 
strong non-ideality of diffusion and viscosity. The results agree with the age old 
wisdom that structure making interactions between the two constituents (here A 
and B) lead to a slower relaxation. The opposite has also been observed for model 
II which has the structure breaking interactions. 

(3) What is also remarkable is that the non-ideality manifests itself in nearly 
opposite ways for the two models. This has been observed for all the transport 
properties and hence reflected in the re-entrance behavior also. 

(4) The theoretical calculations reveal that the main reason for the anomalous 
composition dependence of viscosity lies in the variation of the mean square stress 
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fluctuation (MSSF) with the composition of the mixture. Similarly, for friction, it is 
the Einstein frequency which shows non-monotonic behavior. It is thus fair to say 
that the anomalies have both a structural and a dynamic origin. 

This work was supported in parts by the Department of Science and Technology, 
India and CSIR, New Delhi, India. G. Srinivas thanks CSIR for a research fellowship. 
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Figure Captions 



Figure 1 The simulated values showing the re-entrant behavior of the relaxation 
times Tj, are plotted against simulated viscosity for model I. Filled circles represent 
ta while the open circles tb- The direction of the arrow shows the increasing solute 
(B) composition in both the cases. T* — 1.0, p* — 0.85. 

Figure 2 The diffusion coefficients obtained from MD simulation and mode 
coupling theory are plotted for model 1. Filled and open dries represent the solvent 
and solute diffusion coefficients obtinaed from simulations, respectively. Full and 
dashed lines show the MCT results. T* = 1.0, p* = 0.85. 

Figure 3 The diffusion coefficients obtained from MD simulation and mode 
coupling theory are plotted for model II. Filled and open cirles represent the solvent 
and solute diffusion coefficients obtained from simulations, respectively. Full and 
dashed fines show the MCT results. T* = 1.24, p* = 0.85. 

Figure 4 The composition dependence of viscosity obtained from MD simula- 
tions (symbols) and mode coupling theory (lines) for both the models. Filled (open) 
circles give simulation results for model I (model II). The fines give the theories. 
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